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HOT JALISATION APPARATUS AND METHODS 
FIELD OF INVENTION 

This invention relates to apparatus, methods and computer program code 
for.equaUsation in communication systems in which a receiver receives signals from a 
plurality of transmit antennas, and in particular where the pluraUty of transmit antennas 
is associated with a single transmitter. 

BACKGROUND OF THE INVENTION 

A particular problem arises in a communications links where a transmitter with more 
than one transmit antenna is employed since signals received from different transmit 
antennas interfere with one another. This results in so-called multi-stream interference 
(MSI) and causes decoding difficulties. The potential advantage, however, is greatly 
increased throughput (that is, a higher bit rate) for such a communications link, hi this 
type of MIMO (Multiple-input MuUiple-output) communication link the "input" (to a 
matrix channel) is provided by the transmitter's plurality of transmit antennas and the 
"output" (from a matrix channel) is provided by a plurahty of receive antennas. Thus 
each receive antenna receives a combination of signals from all the transmitter's 
transmit antennas which must be unscrambled. 

A typical wireless network comprises a plurality of mobile terminals (MT) each in radio 
communication with an access point (AP) or base station of the network. The access 
points are also in communication with a central controller (CC) which in turn may have 
a link to other networks, for example a fixed Ethemet-type network. Until recently 
considerable effort was put into designing systems so as to mitigate for the perceived 
detrimental effects of multipath propagation, especially prevalent in wireless LAN 
(local area network) and other mobile communications environments. However the 
described work G. J. Foschini and M.J. Gans, "On limits of wireless communications in 



a fading environment when using multiple antennas" Wireless Personal 
Communications vol. 6, no.3, pp.31 1-335, 1998 has shown that by utiUsing multiple 
antenna architectures at both the transmitter and receiver (a so-called multiple-input 
multiple-output (MIMO) architecture) vastly increased channel capacities are possible. 
Attention has also turned to the adoption of space-time coding techniques (in OFDM, 
space-frequency coding) for wideband channels. Typically channel state information 
(CSI) for maximum likelihood detection of such coding is acquired via training 
sequences and the resulting CSI estimates are then fed to a Viterbi decoder. 

A techiiique fdr space-time code detection in a MIMO system based upon the use of 
periodic pilot sequences and interpolation filters is described in A. Naguib, V. Tarokh, 
N Seshadri and A. Calderbank "A space-time coding based model for high data rate 
wireless communications" IEEE J-SAC vol. 16, pp. 1459-1478. Oct 1998. 

Figure 1 shows a MIMO communication system 100 in the context of which the 
technique described in Naguib et al operates. An information source 101 provides an 
information symbol s(l) at time 1 to a space-time encoder 102 which encodes the symbol 
as N code symbols c,(l) C2(l). • .., cnO), each of which is transmitted simultaneously from 
one of transmit antennas 104. A plurality M of receive antennas 106 receives 
respectively signals ri(l), . ...tm(1) which are input to receiver 108. The receiver 108 
provides on output 1 10 an estimate s(l) of the encoded transmitted symbol § (1). There 
is a plurality of channels between the transmit and receive antennas, for example all 
channels with two fransmit antennas and two receive antennas. The technique described 
in Naguib et al uses periodic pilot sequences in the transmitted signal to estimate the 
time varying responses of these channels. 

Third generation mobile phone networks use CDMA (Code Division Multiple Access) 
spread spectrum signals for communicating across the radio interface between a mobile 
station and a base station. These 3G networks are encompassed by the International 
Mobile Telecommunications IMT-2000 standard (wwwjtuint). Collectively the radio 
access portion of a 3G network is referred to as UTRAN (Universal Terrestrial Radio 
Access Network) and a network comprising UTRAN access networks is known as a 
UMTS (Universal Mobile Telecommunications System) network. The UMTS system is 



the subject of standards produced by the Third Generation Partnership Project (3GPP, 
3GPP2), technical specifications for which can be found at www.3gpp.org. Fourth 
generation mobile phone networks, although not yet defined, may employ MIMO-based 
techniques. 

In practical data communication systems multipath within a channel results in 
intersymbol interference (ISI), which is often corrected with a combination of 
equalisation and forward error coding. For example a linear equaliser effectively 
convolves the received data with an inverse of the channel impulse response to produce 
data estimates with ISI substantially removed. An optimal equaliser may employ 
maximum likelihood (ML) sequence estimation or maximum a priori estimation (MAP), 
for example using a Viterbi algorithm. Where data has been protected with a 
convolutional code a soft input Viterbi decoder may be employed, usually together with 
data interleaving to reduce the effects of burst errors. Such approaches provide optimal 
equalisation but become impractical as the symbol alphabet size and sequence length (or 
equivalently channel impulse response length) increases. 

Turbo equalisation achieves results which are close to optimal but with substantially 
reduced complexity compared to non-iterative joint channel equalisation and decoding. 
Broadly speaking turbo equalisation refers to an iterative process in which soft 
(likelihood) information is exchanged between an equaliser and a decoder until a 
consensus is reached. The effect of the channel response on the data symbols is treated 
similarly to an error correction code and typically a soft output Viterbi algorithm 
(SOVA) is used for both. Again, however, such techniques are impractically complex 
for large delay spreads and symbol alphabets, particularly as several processing 
iterations may be needed to achieve convergence for a single data block. These 
difficulties are significantly exacerbated where signals firora more than one transmit 
antenna must be disentangled and equalised, with a different channel response for each 
transmit antenna or transmit-receive antenna pair. 

A paper by Tuchler et al. (Minimum Mean Squared Error Equalization Using A-priori 
Information, Michael Tuchler, Andrew Singer, Ralf Koetter, IEEE Transactions on 
Signal Processing, vol. 50, pp. 673 - 683, March 2002) describes a simplified approach 
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to turbo-equalisation where a single transmit antenna is employed. In this paper the 
conventional MAP equaliser is replaced by a linear equaliser (that is by a linear or 
transversal filter) with filter coefficients which are updated using a minimum mean 
square error (MMSE) criterion evaluated over both the distribution of noise and the 
distribution of symbols. A linear estimate x„ of a transmitted symbol x„ is determined 
using an observation z„ via the equation x„ = a„" z„ + b„ where superscript " denotes the 
Hermitian operator and a„ and b„ are the coefficients of the estimator (strictly the 
estimate should be termed afifine rather than linear because of the constant b„). The 
coefficients are chosen to minimise the MSE cost £(|x„ - x„ f ) where E( ) denotes a 
mean or expectation value. As information is fed back to the equaUser fi-om the error 
correction decoder the filter coefficients change with time and are thus recomputed for 
each data symbol to be estimated. A related technique is described in WO 02/15459. 

The contents of the Tuchler et al. paper are helpfiil for understanding the present 
invention, which builds upon and extends this work and, in particular, sections II and III 
of this paper are specifically hereby incorporated by reference. 

Another paper, Tetsushi Abe and Tad Matsumoto, "Space-Time Turbo Equalization and 
Symbol Detection in Frequency Selective MIMO Channels" in: Proc. Veh. Techn. 
Conference, IEEE VTS 5*. Vol. 2. pg 1230-1234, 2001 describes the appUcation of 
turbo equaUsation to a MIMO system with a plurality of users (transmitters). However 
the simplified approach described in this paper is only suitable for BPSK and not for 
other (higher) modulation schemes used, for example, by wireless LAN networks. 
Therefore improved techniques, and in particular improved algorithms for systems with 
multiple antenna transmitters are still desirable. 

Other background prior art may be found in US 2002/01 10188 which also employs 
adaptive linear filtering in SISO (soft-in soft-out) turbo equalisation, in Iterative (Turbo) 
Soft Interference Cancellation and Decoding for Coded CDMA, Xiaodong Wang, H. 
Vincent Poor, IEEE Trans. Comms. Vol 47, No. 7, July 1999, pp. 1046-1061, and m 
Iterative Receivers for Multiuser Space-Time Coding Systems, Ben Lu and Xiaodong 
Wang, IEEE Journal On Selected Areas in Communications, Vol. 18 No. 1 1, Nov 2000. 



5 

pp. 2322-2335. In these documents the references to "soft" information are to 
information relating to the likelihood of a particular bit, symbol or signal value as 
opposed to, for example, hard information resulting from a hard bit decision defining a 
bit as a logic one or zero. More complex trellis-based and Bayesian-based SISO 
equaliser components have also been utilized for iterative turbo-equalization where the 
equalizer and decoder iteratively exchange transmitted symbol likelihood information to 
improve the bit-error-rate of the receiver. 

Filter-based equaUsers using MMSE or zero-forcing criteria, Bayesian-based equalisers 
and trelUs-based equalisers such as (MAP) and Maximum Likelihood Sequence 
Estimation (MLSE) equalisers have all been employed for wideband MIMO systems to 
mitigate both ISI and MSI. It is also known to use iterative turbo equalization 
technology for equalisation in MIMO systems. However turbo equaUsers employing 
treUis or Bayesian-based SISO equalisation have a high computational complexity, and 
this increases exponentially with the length of channel response and the number of 
possible states/levels of the modulation employed. The simpUfied turbo equalisation 
schemes described in Tuchler et al., WO 02/15459 and US 2002/01 10188 consider 
systems with a plurality of receive antennas but are incapable of equalisation in systems 
with more than one transmit antenna. It will be appreciated that with multiple receive 
antennas and a single transmit antenna each receive antenna is associated with a single 
channel whereas the situation is much more complex with multiple transmit antennas 
since a combination of signals from different transmit antennas is received. 

There is therefore a need for a reduced computational complexity turbo equahzation 
scheme suitable for mitigating both MSI and ISI in a MIMO communications system. 



STTMMARY OF THE TNVENTION 

According to a first aspect of the invention there is therefore provided a soft-in-soft-out 
(SISO) equaliser for use in a receiver of a communications system employing a plurality 
of transmit antennas, the equaliser comprising: at least one received signal input for 
inputting a received signal; a plurality of likelihood value inputs, one for each transmit 
antenna, for inputting a plurality of decoded signal Ukehhood values from a SISO 
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decoder; a processor configured to determine fi-om said plurality of signal likelihood 
values an estimated mean and covariance value for a signal fi-om each of said transmit 
antennas; an expected signal determiner coupled to said processor to determine an 
expected received signal value using said mean values; a subtractor coupled to said 
received signal input to subtract said expected received signal value from said received 
signal to provide a compensated signal; a filter coupled to said subtractor to filter said 
compensated signal to provide a plurality of estimated transmitted signal values, one for 
each said transmit antenna; a filter coefficient determiner coupled to said processor to 
determine coefficients of said filter using said covariance values; and an output stage 
coupled to said filter to output a plurality of transmitted signal likelihood values, one for 
each said transmit antenna, for input to said SISO decoder. 

The equaliser may be used in a receiver of a MIMO communications system in which 
case a plurality of received signal inputs is provided and a MIMO filter is employed. 
Broadly speaking use of a multiple output or MIMO linear (affine) or transversal filter 
provides a computational complexity which in embodiments increases only Unearly 
with increasing channel impulse response (CIR) length and size of signal constellation. 
With the aid of soft information from the decoder both ISI fix>m the dispersive channels 
and MSI from the signals from the multiple transmit antennas may be jointly mitigated 
(in embodiments of the invention some residual ISI may be present at the equaUser 
outputs but MSI may be substantially removed). Embodiments of the invention thus are 
able to provide a low computational complexity iterative turbo equalisation solution for 
MIMO systems where conventional techniques are impractical. Embodiments of tiie 
equaliser are particularly usefiJ for relatively highly dispersive channels. Any suitable 
conventional decoder may be employed with the equaUser including, but not limited to, 
convolutional, block, tiiibo, low density parity check or space-time tireUis decoders, or a 
serially or parallel concatenation of these, depending upon the encoder used at the 
transmitter. 

In embodiments the expected received signal value comprises an ISI (and/or MSI) 
component of the received signal, for example derived by applying a channel response 
(H) to an expected or estimated transmitted signal (£(X)). Preferably the expected 
received signal value is adjusted by an amount dependent upon expected values of 



estimated transmitted signal values derived from the signal likelihood values from the 
SISO decoder, for example expected symbol interference values. For an n-th 
transmitted symbol interference to the signal x„' from the i-th transmit antenna may be 
determined as described later from e'£(xj) where e' is the i-th column of a symbol 
interference matrix S. 

Preferably the filter coefficients are determined according to a mean square error cost 
fimction, in particular a MMSE function and a filter gain controller is provided to 
confrol the gain of the filter dependent upon co variance values of estimated transmitted 
signal values, as derived from this signal likelihood values. Agaiii the mathematical 
basis of this is described later. 

To fiirther reduce the computational complexity the filter coefficients may be held 
constant over a packet or block of received symbols rather than recomputed each 
symbol. Similarly the filter coefficients may be computed m the frequency rather than 
in the time domain, for example in a single-carrier modulation system by transforming 
a packet or block of received symbols to the frequency domain for processing in parallel 
and then afterwards transforming back into the time domain. Alternatively, in an 
OFDM (Orthogonal Frequency Division Multiplexed) system, where the symbols are in 
the frequency domain, a transformation such as an IFFT (inverse Fast Fourier 
Transform) is not necessary after frequency-domain equahsation. 

A receiver may be provided incorporating both the equaliser and the decoder and 
preferably, to improve the mdependence of the transmitted data, also a deinterleaver 
following the equaliser and an interleaver to interleave the likelihood value inputs to the 
equaUser from the decoder. Such a receiver may be incorporated in, for example, a 
mobile terminal, base station or access point of a communications system witii more 
than one transmit antenna, for example for increased data rates and/or, where transmit 
diversity is employed, for increased robustness. 

According to another aspect of the invention there is provided a method of equaHsing 
received data in a communications system employing a receiver with at least one 
receive antenna and a transmitter with a pluraUty of transmit antennas, the method 
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comprising: inputting at least one received signal from said at least one receive antenna; 
inputting a plurality of decoded signal likelihood values, one for each transmit antenna, 
from a SISO decoder; determining from said plurality of signal likelihood values an 
estimated mean and covariance value for a signal from each of said transmit antennas; 
determining an expected received signal value using said mean values; subtracting said 
expected received signal value from said received signal to provide a compensated 
signal; determining coefficients for a filter using said covariance values; filtering said 
compensated signal in accordance with said coefficients to provide a plurality of 
estimated fransmitted signal values, one for each said transmit antenna; and outputting a 
plurality of transmitted signal likelihood values, one for each said transmit antenna, 
derived from said estimated transmitted signal values, for input to said SISO decoder. 

According to a fiirther aspect of the invention there is provided a SISO MIMO turbo- 
equaliser configured for use with a SISO decoder to equaUse data from a multiple 
antenna transmitter, the equaUser comprising a multi-dimensional transversal filter 
having a plurality of soft inputs and providing a plurality of soft outputs, the equaliser 
being configured to receive soft information firom said SISO decoder and to use said 
soft information to adjust coefficients of said transversal filter to mitigate MSI and ISI. 

hi another aspect the invention provides a method of turbo-equalising data in a receiver 
receiving data from a transmitter configured to transmit data from a plurality of transmit 
antennas simultaneously, the receiver including a SISO decoder, the method 
comprising: inputting a received signal vector comprising a block of received signal 
data including data received from said plurality of transmit antennas; receiving soft 
decoded information from said decoder; processing said received data using a linear 
filter having a plurality of soft inputs derived from said received signal vector to 
provide a plurality of soft equaliser outputs, one for each transmit antenna, to said SISO 
decoder; and adjusting coefficients of said linear filter in response to said soft decoded 
information to mitigate MSI and ISI. 



The invention also provides a method of equalising data in a receiver of a 
communications system with a plurality n, of transmit antennas the method comprising: 



inputting a received signal vector Z„ comprising a block of received signal data at an 
index «; inputting from a SISO decoder a set of bit likelihood values l(c'„j )t i = 1 to m, 
j = 1 to m where c\ j, denotes a portion of an m-bit codeword symbol at an index «; 

determining expectation e(x^„ ) and covariance cov(x; ) values for estimated 
transmitted signal values x[ using the values l[c[ j\ where x[ denotes an estimated 
value of a signal transmitted from a transmit anterma / at an index «; determining an 
updated estimated transmitted signal value x[ at index n for each transmit anteima / 

using =A::ff[z„ -(£(Z„)-^e'£:(x:))] where E{Z„)=1AE{X„) and H is a channel 
impulse response matrix and E(K,J is an expectation value of an estimated transmitted 
signal vector X„ at index n derived from values e[x[ ), where e' is the ith column of a 
matrix S and S = h[o„,,(^,,,_,)„, 0„^,^,„, Y where matrix I,,, is an / x identity 

matrix and 0. is a matrix in which each element is substantially zero, is the ith 
column of a filter matrix F„ and F„ is derived fi-om cov(x„ , x J , and where is a 
scalar filter gain; and outputting equalised likelihood values (< derived from 
values x\ . 

The transmitted symbol index n may index time and/or firequency, depending upon the 
modulation scheme employed. 

The skilled person will recognise that the above-described equahsers and methods may 
be implemented using and/or embodied in processor control code. Thus in a further 
aspect the invention provides such code, for example on a carrier mediimi such as a 
disk, CD- or DVD-ROM, programmed memory such as read-only memory (Firmware) 
or on a data carrier such as an optical or electrical signal carrier. Embodiments of the 
invention may be implemented on a DSP (Digital Signal Processor), ASIC (Application 
Specific Integrated Circuit) or FPGA (Field Programmable Gate Array). Thus the code 
may comprise conventional program code, or micro-code, or, for example, code for 
setting up or controlling an ASIC or FPGA. In some embodiments the code may 
comprise code for a hardware description language such as Verilog (Trade Mark) or 
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VHDL (Very high speed integrated circuit Hardware Description Language). As the 
skilled person will appreciate, processor control code for embodiments of the invention 
may be distributed between a plurality of coupled components in communication with 
one another. 

BRIEF DESCRIPTION OF THR DRAWINGS 

These and other aspects of the invention will now be further described, by way of 
example only, with reference to the accompanying figures in which: 

Figure 1 shows a known MIMO space-time coded communications system; 
Figure 2 shows a block diagram of a MIMO channel model; 

Figure 3 shows a flow diagram of a MMSE MIMO turbo equalisation procedure 
according to an embodiment of the present invention; 

Figure 4 shows a block diagram of a soft-in-soft-out filter-based MIMO equalizer 
according to an embodiment of the present invention; 

Figure 5 shows a MIMO communications system including a receiver employing the 
MIMO equaliser of Figure 4; 

Figure 6 shows a block diagram of a fi-equency domain soft-in-soft-out filter-based 
MIMO equalizer according to a fiarther embodiment of the present invention; 

Figure 7 shows a receiver incorporating an equaliser configured to operate in 
accordance with an embodiment of an aspect of the present invention; and 

Figure 8 shows bit error rate performance against received signal to noise ratio for a 
turbo-equaUzer and according to an embodiment of the present invention, for varying 
numbers of equalization-decode iterations. 
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nFTATLED PFSCRIPTION 



Referring to Figure 2, this shows a block diagram 200 of a MIMO channel model. A 
multi-stream transmitter has first plurality n, of transmit antennas 202 and transmits 
respective symbols x„' , x„^. . at time n which comprise "inputs" to a matrix 
channel 206. A plurality Ho of receive antennas 204 provides "outputs" firom the 

n, X rio MIMO matrix channel in the form of received signals Zn\-z^ z"„° • There 

is a channel with an associated channel response hn'"* between the ith transmit antenna 
and jth receive antenna. The received signal at each receive antenna also includes a 
noise component Wri" . 

In an X no MIMO system the MIMO equalizer aims to provide an estimate of n, 
transmitted data symbols at every signalling instant. Data firom the plurality of transmit 
antennas are transmitted at the same or overiapping times and using at the same or 
overiapping fi-equencies and thus MSI is introduced, as well as ISI firom the dispersive 
wideband channel. In some systems the same or related data is transmitted firom 
different transmit antennas in order to provide redundancy or diversity, hi other 
systems different data streams are transmitted from each transmit antenna, for example 
to provide higher overall data rates. 

Broadly speaking, in embodiments of the invention a turbo equalizer is m the form of a 
multi-dimensional transversal filter. Coefficients of the filter are adjusted according to 
a Minimum Mean Square Error (MMSE) criterion to mitigate the effects of both Inter- 
Symbol Interference (ISI) and Multi-Stream hiterference (MSI) and provide a soft 
output, with the aid of soft information relating to multiple transmitted data signals 
received from a decoder. The equalizer processes signals from all the tiransmit antennas 
in parallel (although potentially in a time-multiplexed manner in processor-based 
embodiments) is thus able to "detect" and equaUse amultistiream signal. The soft or 
likelihood information from the decoder is used to determine tiie covariance and mean 
of tiie multistream transmitted signal for evaluating the coefficients of the equatizer and 
also to determine a mean or expected value of the received signal for MMSE based 
detection. 
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In embodiments the SISO equalization may be perfomed in either or both of the time 
and frequency domains. Initially a time-domain implementation will be described; 
details of an alternative frequency-domain implementation will be given later. 

The elements of the equaliser will generally be implemented by a digital signal 
processor and the structure and operation of a first, time-domain embodiment of the 
equaliser will therefore first be described in mathematical terms. 

Consider first a signal' transiiiission ihbdel for a linear, dispersive and noisy 
communication system with n, inputs and /i^ outputs (that is a x MIMO 
system). The complex-valued baseband equivalent signal model can be represented as 
by equation (1) below: 



A:=0 



(1) 



where L is the channel length in symbol periods, where x„=[z\ zl ■■• zl°] is the 
signal received by the «o receive antennas, where w „ = [w^ wl ■■■ w^" f is the 



noise 



component, where x„=k • •• x„"' f is signal transmitted by the 

n, transmit antennas, and where n labels a symbol (ie. in this case corresponds to time 
in symbol periods). The dispersive channel matrix is denoted by: 









... hr' 






hi' 


... 






hi'"" 


... hi""" 



Here, h);' represents the ifcth channel tap for a channel link between an rth transmit 
antenna and a yth receive antenna as shown in Figure 2. 
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Consider an embodiment of a transversal filter based equaliser with (integer) N 
coefficients or taps where iy^ = M + A^2 + 1 (M and iV2 in effect specifying the length 
of the causal and non-causal portions of the filter). The filter employs an observation of 
a block of received signal of length N in symbol periods, which is denoted as: 



Z„=HX„+W„ 



(2) 



where 



4- 



n+Nl 



(3) 



is the Ntia X 1 received signal vector, 



(4) 



is the Nho x 1 received noise vector, 



(5) 



is the + L - 1)«/ X 1 transmit signal vector which convolves with the 
Ml x{N + L- 1)«, channel impulse response (CIR) matrix H given by: 



H = 



X-2 



*»0 ^no>^, 



0 



L-2 



0 h, , h 



L-2 



(6) 



Determination of the filter coefficients and of a soft transmitted symbol estimate by 
equaliser will next be considered. 
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An optimum linear estimate of the transmitted symbols x„ = [x^ ■■■ x"' f upon 

receiving Z„ , based on the second-order statistics (that is, means and covariances) has 
the form given in equation (7) below (see, for example, "An Introduction to Signal 
Detection and Estimation", H. Vincent Poor, Second Edition, Springer, ISBNO-387- 
94173-8 at pages 221 to 233): 

i„=^(xj4.A/[Z„-^(Z„)] (7) 

where, as before, £{■) denotes a mean or expectation value and A„ is the set of linear 
filter/equalizer coefficients at symbol (time) n. It can be seen from equation (7) that the 
estimate of the transmitted symbol x„ is the mean offset of the filtered difference 
between the received signal vector Z„ and the ISI and MSI component E{Z„ ) at time n . 

Minimizing the mean square error cost ftmction EQxn - K P)» ^ optimum set of filter 

coefficients A„ is given by: 

A, =cov(Z„,Z„)-'cov(Z„,x„) (8) 

where the covariance operator cov(x,y) = ^(xy") - Eix)E(y^). Equation (7) can be 
expressed in the form: 

x„ = Mx„) + cov(x„,x„)F/[Z„ -E{Z„)] (9) 

whCTe the expectation value of the received signal observed at the equalizer taps, in 
terms of the transmitted signal or an estimate of this signal, is given by: 



E(Z„) = IIE{X„) 



(9a) 
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As will be seen later, E(Z„) may be detennined from the soft decoder output. In 
equation (9) the Miq x «,. fiUer matrix F„ is defined as 

= [R^+HR;^H^rS 
where the covariances are defined as 

Rzz=cov(Z„,Z„), 
^xx= cov(X„,X„), 
R^=cov(W„,Wj. 

The Nrip x rii CIR matrix S introducing the symbol interference firom the n -th 
transmitted symbols is given by: 

S = H[0„^^(jV2+i-l)n, In,x„, 0„^>,;^,„^] (11) 

In equation (1 1) matrix I,,, is an z x i identity matrix and matrix 0,,^. contains all zeros. 

It can be understood from the form of the equalizer (filter) coefficient matrix F„ 
expressed in equation (10), which includes the channel matrix H as defined in equation 
(6), that the equalizer not only performs equalization on the ISI but also mitigates MSI. 

The mean and covariance of the transmitted signal, E{x„) andcov(x„,xJ . which are 
used to compute the filter coefficients and the estimated transmitted signal, may be 
obtained using equations (16) and (17) below. The block-diagonal covariance matrix 
comprises (N+L-l) block-diagonal blocks of cov(x„,xJ as follows: 



16 



GOV (x„_/^2-£+l » ^n-Af 2-1+1 ) 



R XX- 



The cm. matrix H may be obtained from a MIMO channel estimation block in the 
receiver in a conventional manner. Thus in one embodiment known pilot or training 
sequences are periodically inserted into the signal from each transmit antenna and at the 
receiver these known sequences are encoded and provided to a channel estimator 
together with one (or more) input signal streams from the one (or more) receive 
antennas. Providing the training sequences in the signal streams have a low cross- 
correlation coefficient the channel estimator is then able to determine the CIR matrix H 
by, in effect, correlating the encoded known data sequence for each stream with the 
actual received data. An example of such a channel estimator is described in Ye 
Geoffrey Li, "Simplified channel estimation for OFDM systems with multiple transmit 
antennas", IEEE Transactions on Wireless Communications, Vol. 1, No. 1, pg. 67, Jan 
2002, which is hereby incorporated by reference. 

Advantageously the channel estimation may also be performed iteratively by using the 
estimated transmitted symbols Jc^ as additional training symbols. For background 
information relating to this reference may be made to Tetsushi Abe and Tad 
Matsumoto, "Space-Time Turbo Equalization and Symbol Detection in Frequency 
Selective MIMO Channels" in: Proc. Veh. Techn. Conference, IEEE VTS 5*. Vol. 2. pg 
1230-1234, 2001, also hereby incorporated by reference. Thus, for example, multiple, 
deinterleaved outputs from the decoder (one for each fransmit antenna) may be re- 
encoded and provided to the channel estimator similarly to a known fraining sequence. 
In this way, an updated estimated value of H may be obtained on every turbo 
equalization iteration and this may in tum be used to compute the equalizer coefficients. 

For good convergence behaviour the soft transmitted symbol information x'„ provided 
by the equalizer should be independent from the soft fransmitted symbol information 
from the decoder in the iterative structure. At the fransmitter a sequence of m bits, 
where m is the number of bits per symbol, is encoded as a codeword 
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c' .e{+l,-l},7 = l," -,'« for transmit antenna I. The decoder provides soft log 
likelihood ratio information l{c'„ j \j = \,--,m where 

l(c' .) = ln|p(c; ^ = +l) / P{c[ J = -l)} which, from the point of view of the equaliser, 
may be considered as a priori information. To provide good convergence behaviour the 
likelihood information for bits c'„j ,j = l,-,m (corresponding to symbol x'„ ) which are 
available to the equalizer is set to l{cIj f^''^ = 0, y = 1, • • • , m to derive an equation for 
jc^ (During the turbo equalisation these values may also be initiahsed to zero). Then, 
eliminating the dependence of the equalisation on L[c'„j\j = \,-,m provided by the 
decoder, and setting the expectation e[x[ f""'""" =0 and the covariance 
cov(x; , x'„ )^""'""' = 51- , where Sl is an initial value, for example unity (depending 
upon the relative gain of the transmit antennas), for the signal variance of transmitted 
symbols from transmit antenna i , the i th component of the estimated transmitted 
symbol x„ in equation (9) can be expressed as scaled version of F„ (using the matrix 
mversion lemma), as shown in equation (12) below: 



where the filter gain K[ is given by: 

• i^:=(l + fe-cov(<,<)>re'r (13) 

In equations (12) and (13) is the ith column of the filter matrix F„ defined in 
equation (10) and e' is the / th column of the symbol interference matrix S given in 
equation (11). 

Next a reduced complexity method of determining R^z'' is described, for computing 
the fiher matrix F„ via equation (10). In order to avoid performing a direct matrix 
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inversion of the Nno^Nrio covariance matrix Rzz to compute the filter matrix F„ for 
every estimate of symbols x„ , a time-recursive update algorithm may be employed. 
The algorithm is similar to one described in the Tuchler et al. paper (ibid) and exploits 
the existence of common submatrices within a partitioned covariance matrix. The 
matrix Rzz'* at a time n + 1 may be determined by the manipulation of submatrices as 
set out in equation (14): 



R^-'(« + l) = 



— . H 



(14) 



These submatrices may in turn be determined based upon a previous R^^ ^ at a time n , 
as shown by equations (15) below: 



Up u 



Up Up J 



—1 M H 

=[Rj^+HR^H"] 



0, 



{N-\)noyiio 



(15) 



The sizes of the matrices w ^ , cr^^ , , , are (iV - ^n^ x . «o «o . «o ^ «o > 
(JV - 1)mo X "o anti ~ l)"o X - l)"o ' respectively. It can therefore be seen that 
only the matrix inversion of «o x matrices Upand(cr^-W/^"cT^)is required to 
update the matrix R^z'' . The term R^ is the noise covariance and a value for this 
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may be determined, for example, theoretically from the receiver bandwidth (particularly 
where the front end filter has a relatively sharp cut-off; see also 3GPP TS25.215 v5.2.1 
for background on received power), or by a measurement of the level of noise (and/or 
interference) at the receiver, or by a combination of both these techniques. Where the 
noise levels at the receive antennas are similar may approximate to a fraction of I. 

The soft (Ukelihood) value inputs to the equaliser from the decoder will now be 
discussed. 

The soft information provided by the decoder in terms of the likelihood values are 
utilized to provide the mean e{xI ) and covariance cov(xl,xl )of the transmitted 
symbols, which are required to compute in Equation (9). The mean and covariance 
values are obtained as shown below: 



and 



e{<)=11^A<-^>^) (16) 



At the transmitter a sequence of /w bits c is mapped to a symbol x selected from a 
0 = 2" -ary symbol alphabet S = )a„a2.",ae} in which each symbol corresponds 
to an m bit pattern s, = 5,., • • • J, s,, 6 {+ l,-l} . The probabilities of the 
transmitted symbols are given by: 



= n;,, 1 / 2.^ + 5,,,. tanh(z(<,. )/ 2)} 



(18) 
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The soft (likelihood) value outputs from the equaliser to the decoder will next be 
considered. 



At the output of the equalizer, it is assiuned the probability density functions of the 
estimated symbols p{xi 1 = ) are Gaussian with mean /i^^ and variance cT„^,t' . The 
log-Ukelihood ratios of the bits contained in the transmitted symbols are then given by: 



-\E{xl-Mi,,}' 



-In 



-\e{x'„ - mIj^J 



J) 



tlyk 



(19) 



Based on the derivation obtained for the estimate of the transmitted symbol x'„ in 
equation (12) the mean and variance of the equalizer output may then be written: 

=Kf:"[E(zAx[ =aj-(£(zj-e'£(x; 



H 



(20) 



and 



The logarithmic summation of exponentials used in equation (19) to compute the 
equalizer's soft output, hi exp(l, ) , may advantageously be computed usmg the 
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Jacobian logarithmic relationship: 



ln[^g exp(^, )j = J{X„ . J{X„_, •)) 
y(A„A^)»max(A,,A,) + /,(|A, -/l^l) 
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Here f^x) is a logarithmic correction function that can be tabulated in a look-up table 
such as that given in Table 1 below. 





^■4(^)■;^r;■■.:■^•v.;./•:^^:^.:::^^^^^ 


x>3.7 


0.0 


3.7 >a:> 2.25 


0.05 


2.25>x>1.5 


0.15 


1.5>;c>1.05 


0.25 


1.05 > a: > 0.7 


0.35 


0.7 >x> 0.43 


0.45 


0.43>x>0.2 


0.55 


0.2 ^jc 


0.65 



Table 1 



In this way the computational complexity of technique may be reduced by replacing the 
use of logarithmic and exponential functions in the likelihood calculation with a simple 
look-up table and "maxQ" function. 



To further reduce the computation complexity, the equalizer coefficients, that is terms 
F and K' defined in equations (10) and (13) respectively, may be set to be non- 

1 ^ 

varying with time n. An average of the covariance matrix R = — X » ^« ) 

and cov(jc^,<)= — ^ cov{x'„ , x[ ) over a received symbol packet of P symbols, may 

instead be used to compute the coefficients in Equation (10) and (13). These time 
invariant equalizer coefficients F and ^' may then be used to facilitate a reduced 
complexity computation of the transmitted symbol x'„ in Equation (12). 

Figure 3 shows a flow diagram of an implementation the above-described turbo 
equalisation procedure. System parameters are determined or chosen at step S800 
including, for example, the symbol alphabet 5 and the number of filter coefficients 
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m+N2(Nl andN2 are typically chosen such that N\=N1=L). 
however, this is not part of the actual equaUsation procedure. 



Strictly speaking. 



At step S802 equaliser parameters are initialised, for example by setting 
co\{xi X) = ^li ' l[c[_j )= 0 (a log likelihood ratio L(.) of 0 implying c=+l and 
c=-l are equiprobable), then at step S804 a block of received signal Z„ is input from 
one or more receivers and, at step S806, a estimated CIR matrix H is input from a 
channel estimator. 

At step S808 the mean e[x'„ ) and covariance co^x'„,x\, )of the fransmitted symbols are 
obtained, according (after the first iteration) to the a priori information of the 
transmitted bits from the SISO decoder L[c[j),i = = l." -'" (where m is the 

number of bits per symbol) using equations (16), (17) and (18). Then, at step S812, 
using the mean e[x'„) and covariance cov{x[,x[) values of the fransmitted symbols the 
equalizer coefficients F„ and K'„ are obtained (using equations (10) and (13) 
respectively) and the mean of the received signal £(Z„) is determined using equation 
(9a). The estimated transmitted signal x„ is then given by equation (12), and the 
"exfrinsic" information on the fransmitted symbols, , which will provide a soft 

input to the SISO decoder from the equaUser, is obtained using equation (19) given the 
estimated fransmitted symbols and the a priori information relating to the fransmitted 
bits. 

Deinterleaving is performed at step S814, corresponding to the interleaving performed 
after encoding and before transmission of the data at the fransmitter. Thus the extrinsic 
information from the equalizer is deinterleaved before passing to the decoder. Then at 
step 8816 the extrinsic information relating to the transmitted bits is used as the a priori 
knowledge at the decoding stage to provide the decoder's extrinsic information. The 
deinterleaved signals are decoded in accordance with the encoding used at the 
transmitter. A check is then made, at step S818, to determine whether a termination 
criterion has been reached and, if it has not, the procedure loops back to reinterleave the 
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extrinsic information from the decoder and to pass this back to the equalizer via step 
S808. The termination criterion may comprise, for example, a predetermined number of 
iterations or a determination of whether the decoder soft output, that indicates the 
reliability of the decoded data, is more than a threshold value. If the termination 
criterion has been met the equalized and decoded bits at the final iteration are provided 
as an output estimate of the transmitted data bits. 

Figure 4 shows a block diagram of a soft-in-soft-out filter-based MIMO equahzer 400 
configured to operate according to the procedure of Figure 3. 

Referring to Figure 4 the equaliser comprises an input 402 to receive a plurality of 
received signal blocks from a corresponding plurality of receive antennas, these 
received signal inputs providing a first set of inputs to a subtractor 404. A second set of 
inputs to the equaHser 400 comprises a set of soft or likelihood values 406 from a soft- 
in-soft-out decoder, to provide a priori transmitted bit information to the equaUser. 
This information is processed by block 405 to determine E(xi ) and coy(xl,x'„ ) and 
these (estimated) mean (or expectation) and covariance values are used by block 408 to 
perform a calculation (E(Z„)-e'E(x'„)) to determine an expected set of received signal 
values. This provides a second set of inputs 410 to subtractor 404 for subtraction fix>m 
the first set of inputs 402 to provide a "compensated" signal 412 to a Imear (affine) or 
transversal filter 414. The soft likelihood values 406 from SISO decoder are used by a 
filter coefficient calculation block 416 to determine a set of coefficients for filter 414. 
A MIMO channel estimator 41 8 provides a channel estimate input 420 to this filter 
coefficient calculation block 416. The output of filter 414 comprises a set of estimated 
transmitted signal values 422, one for each transmit antenna, which are in tiim provided 
to a soft decision block 424 which provides a corresponding pluraUty of outputs 426 
comprising soft (or likeUhood) transmitted bit values which are suitable for providing a 
soft input to the SISO decoder. 

Figure 5 shows a MIMO communication system 500 including a MIMO tiransmitter 502 
and a MIMO receiver 504 communicating via a MIMO channel 506. The receiver 504 
incorporates a MIMO SISO (MMSE) equahser 400, for example as described with 
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reference to Figure 4. The transmitter 502 has a data input 508 providing an input to a 
space-time and/or channel encoder 510 (in other embodiments other coding schemes 
may be employed). The encoder 510 provides a plurality of outputs 512 to an 
interleaver 514 which, in turn, provides signals to a plurality of rf output stages and 
thence to a corresponding plurality of transmit antennas. In Figure 5, for clarity these rf 
stages and transmit antennas are not shown explicitly but instead a plurality of outputs 
516 comprising the plurahty of transmitted signals is shown providing an output from 
interleaver 514 which, via MEMO channel 506, provides a (different) plurality of inputs 
to receiver 504. It will be recognised that the communication system of Figiu-e 5 
employs a single transmitter to provide a plurality of transmit output streams, either for 
redundancy or increased bit rate. The transmissions from the plurality of transmit 
antennas may, for example, share the frequency or overlap in frequency and/or overlap 
in time. This is different to a communication system employing a plurality of users 
with frequency and/or time domain controlled access in which, generally speaking, it is 
preferred to assign different frequencies and/or time slots to different uses. 

Continuing to refer to Figure 5, a plurality of receive antennas coupled to a 
corresponding plurality of rf receiver front ends (not shown for clarity in Figure 5) 
provides a plurahty of inputs to MIMO SISO MMSE equaliser 400. The soft output 
426 from the equaliser 400 are deinterleaved by a deinterleaver 518 and then provided 
to a space-time/channel decoder 520. The decoder accepts a plurality of inputs, one for 
each signal stream from a transmit antenna, and provides a corresponding plurality of 
outputs 522a, b which are either provided to a bit interleaver 524 and retumed to 
equaUser 400 for a fiirther equalisation-decoding iteration or, if a termination criterion 
has been reached, output as estimated data. It will be appreciated that equaliser needs 
complete transmitted symbols from decoder 520, that is where, for example, error check 
bit such as parity bits have been included transmit symbols including these parity bits 
should be provided to the equaliser. The MIMO channel estimator 418 of Figure 4 (not 
shown in Figure 5) may receive a set of inputs from bit interleaver 524 and output an 
estimate of H to equahser 400. 

Referring now to Figure 6 this shows a block diagram of a soft-in-soft-out filter-based 
MIMO equalizer 600 configured to operate in the frequency domain. The main 
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elements of equaliser 600 are similar to those of Figure 4 and are indicated by like 
reference numerals. The main differences are that calculation blocks 408 and 416 
perform frequency domain calculations and that subtractor 404 and filter 414 operate in 
the firequency rather than the time domain (though, for simpUcity of comparison, the 
same reference numerals have been employed and the channel estimator has been 
omitted). Thus the main additions are a serial-to-parallel conversion block 602 and fast 
Fourier transform (FFT) block 604 prior to subtractor 404 and filter 414, and an inverse 
fast Fourier transform (FFT'^) 606 and a parallel-to-serial conversion block 608 
following filter 414- These blocks convert the received signals firom the time to 
firequency domain for processing and back to the time domain after processing. 

Background information relating to SISO equalisation in the fi-equency domain can be 
found in Tuchler, M., Hagenauer, J.,'^Turbo equalization' using firequency domain 
equalizers" in: Proc. Allerton Conference, Monticello, IL, USA, Oct 2000, to which 
reference may be made. It will be appreciated, however, that this reference only 
discloses fi-equency domain turbo equaUsation for a single-transmit-antenna-single- 
receive-antenna system. Other background information can be foimd in M. V. Clark, 
"Adaptive Frequency-Domain EquaUzation and Diversity Combining for Broadband 
Wireless Communications", IEEE JSAC, vol. 16, no. 8, pp. 1385-1395, Oct. 1998, and 
D. Falconer, S. L. Ariyavistakul, A. Benyamin-Seeyar, B. Eidson, "Frequency Domain 
Equalization for Single-Carrier Broadband Wireless Systems", IEEE Communications 
Magazine, pp. 58-66, April 2002. 

Continuing to refer to the example of a frequency domain implementation of the 
equaliser shown in Figure 6, the signals from the receive antenna(s) are transformed into 
the frequency domain by Fast Fourier Transform (FFT) on a per FTT block/packet basis 
and the spectra of the received signals are equalized by multiplying the spectrum of 
each branch with the frequency domain coeflScients of the equalizer. Frequency domain 
equalization again reduces the computational complexity of the equaUzation task by 
processing an entire FFT block/packet of received symbols at once in the frequency 
domain instead of symbol-by-symbol in the time domain although it will be appreciated 
that because of this the equalizer coefficients are time-invariant per FFT block. 
Furthermore a cyclic extension (for example, prefix) should be added to the transmitted 
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symbol block/packet before transmission (to provide an integral number of cycles 
within the FFT interval) although this is discarded at the receiver. The frequency 
domain implementation of the turbo equaliser is suitable for both time and/or frequency 
domain coded data, and may be used, for example, with MIMO OFDM transmissions. 

Referring now to Figure 7, this shows an embodiment of a receiver 700 incorporating a 
equaliser configured to operate as described above, in particular to operate as described 
with reference to Figures 3 and 4. 

Receiver 700 comprises one or more receive antennas 702a, b (of which two are shown 
in the illustrated embodiment) each coupled to a respective rf front end 704a, b, for 
example similar to the rf front end of Figure 2a, and thence to a respective analogue-to- 
digital converter 706a,b and to a digital signal processor (DSP) 708. DSP 708 will 
typically include one or more processors 708a (for example, for a parallel 
implementation of filter 414) and some working memory 708b. The DSP 708 has a 
data output 710 and an address, data and control bus 712 to couple the DSP to 
permanent program memory 714 such as flash RAM or ROM. Permanent program 
memory 714 stores code and optionally data structures or data structure definitions for 
DSP 708. 

As illustrated program memory 714 includes SISO equalisation code 714a comprising 
(£(Z„ ) - e'E(xl )) calculation code, subtraction code, filter coefficient calculation code, 
linear (transversal) filter code and soft decision output code to, when running on DSP 
708, implement the corresponding functions as described in detail above. Program 
memory 714 also includes MIMO channel estimation code 714b to provide a MIMO 
cm. estimate H, de-interleaver code 714c, interleaver code 714d, and space-time (or 
other) SISO decoder code 714e. Examples of suitable code 714e are well known to 
those skilled in the art. Optionally the code in permanent program memory 714 may be 
provided on a carrier such as an optical or electrical signal carrier or, as illustrated in 
Figure 7, a floppy disk 716. The data output 710 from DSP 708 is provided to fiuther 
data processing elements of receiver 700 (not shown in Figure 7) as desired. These may 
include a block error decoder such as a Reed-Solomon decoder (although this could be 
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part of the turbo decoder), and a baseband data processor for implementing higher level 

protocols. 

The receiver front-end will generally be implemented in hardware whilst the receiver 
processing will usually be implemented at least partially in software although one or 
more ASICs and/or FPGAs may also be employed. The skilled person will recognise 
that all the functions of the receiver could be performed in hardware and that the exact 
point at which the signal is digitised in a software radio will generally depend upon a 
cost/complexity/power consumption trade-ofif. 

Figure 8 shows bit error rate (BER) performance against received signal to noise ratio 
(El/No) for an embodiment of the tmbo-equalizer. for varying numbers of equalization- 
decode iterations. In figure 8 curve 800 shows a theoretical result for a single transmit 
single receive antemia system with a non dispersive chamiel and the remaining curves 
show the 8PSK performance of an MMSE-based MIMO SISO turbo-equalizer for a 
two-transmit-antenna transmitter, two-receive-antemia receiver MIMO system with a 

Ai-' =0.707,A,'-' = 0.707, =0-707, A,'-' =0.707, 
dispersive channel of =0.707,;i,^' =0J07,K' =0.9579,^f'^ =0.2874' 
TTie system utiUzes half-rate (5.7) convolutional coding and the equalizer parameters are 
Nl=9 mdN2=5. More specifically curves 802. 804, 806. 808 and 810 show, 
respectively, curves for one. two. three, four and five turbo equalisation iterations. It 
can be seen that performance increases significantly with the number of iterations and 
approaches that of a non-dispersive channel, although with large numbers of iterations 
the incremental performance increase is reduced. The reduction in incremental increase 
can be used to determine a termination criterion for a potential application, for example 
limiting to 2, 5 or 10 iterations depending, for example, upon speed-processing power 
trade-offs. 

Applications of the invention have been mainly described in the context of a MIMO 
system with time domain coding but embodiments of the invention are also useful in 
frequency domain coded systems such as MIMO-OFDM (Orthogonal Frequency 
Division Multiplexed) systems. For example the invention may be employed with the 
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European Hiperlan/2 or US IEEE 802.11a standards for 54Mbps wireless networks. 
Embodiments of the invention may also be employed in non-wireless applications such 
as magnetic or optical disk drive read head circuitry where, for example, multiple layers 
of a disk in effect act as multiple transmitters, one or more heads receiving read data 
influenced by "transmitted" signals from more than one layer. 

No doubt many other effective alternatives will occur to the skilled person. It will be 
understood that the invention is not limited to the described embodiments and 
encompasses modifications apparent to those skilled in the art lying within the spirit and 
scope of the claims appended hereto. 



